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Abstract. We study generalised anomalous diffusion processes whose diffusion 
coefficient D{x,t) ^ Do\x\°‘t^ depends on both the position x of the test particle 
and the process time t. This process thus combines the features of scaled Brownian 
motion and heterogeneous diffusion parent processes. We compute the ensemble and 
time averaged mean squared displacements of this generalised diffusion process. The 
scaling exponent of the ensemble averaged mean squared displacement is shown to 
be the product of the critical exponents of the parent processes, and describes both 
subdiffusive and superdiffusive systems. We quantify the amplitude fluctuations of 
the time averaged mean squared displacement as function of the length of the time 
series and the lag time. In particular, we observe a weak ergodicity breaking of 
this generalised diffusion process: even in the long time limit the ensemble and time 
averaged mean squared displacements are strictly disparate. When we start to observe 
this process some time after its initiation we observe distinct features of ageing. We 
derive a universal ageing factor for the time averaged mean squared displacement 
containing all information on the ageing time and the measurement time. External 
confinement is shown to alter the magnitudes and statistics of the ensemble and time 
averaged mean squared displacements. 


PACS numbers: 05.40.-a., 02.50.Ey, 87.I0.Mn 

1. Introduction 

Deviations from Brownian motion are qnite nbiqnitons in a large variety of complex 
systems. Mostly snch anomalons diffnsion processes are characterised by a scaling form 
of the mean sqnared displacement (MSD), 

( 1 ) 

The magnitnde of the anomalons scaling exponent k , distingnishes snbdiffnsion for 
0 < /t < 1 and snperdiffnsion k > 1 [DEIEIII]. Examples for anomalons inclnde the 
relative diffnsion of tracers particles in fnlly tnrbnlent or weakly chaotic [U] systems 
as well as in gronndwater aqnifers [7]. Snbdiffnsion of charge carriers in amorphons 
semiconductors was originally analysed some 40 years ago [8] bnt is now regaining 
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attention in the stndy of polymeric semiconductors [9]. However, the main current 
impetus for the study of anomalous diffusion processes is due to modern spectroscopic 
tools such as fluorescence correlation spectroscopy [10] and single particle tracking of 
submicron particles m- By these methods pronounced deviations from normal diffusion 
(k = 1) were found in structured and crowded liquids |I2l HSl HU US] and in living 
biological cells HU [HI HI]. Similarly, anomalous diffusion is observed in supercomputer 
studies of, inter alia, flexible networks [12] and lipid membranes |2Uj . 

The non-Brownian scaling of the MSD of the tracer particles may originate from 
a range of physical mechanisms. These include the subdiffusive motion on geometric 
fractals such as percolation clusters close to criticality p Hsi ini or the anomalous 
diffusion in Lorentz gases Another important class of anomalous diffusion models 
are continuous time random walks (CTRWs), in which the test particle’s motion is 
interrupted by random waiting times [23] . If the distribution of these waiting times 
is scale free, subdiffusion emerges [8], a phenomenon closely related to quenched trap 
models with exponentially distributed trap depths fflEU. The third class represents 
processes, in which the particle is driven by fractional Gaussian noise which is long 
range temporally correlated: the fractional Brownian motion [2S| or associated fractional 
Langevin equation motion m are connected to viscoelastic environments and represent 
the motion of tagged monomers in a Rouse chain [27] or a tracer particle in a single file 
|28j . A more complete overview of anomalous diffusion models provides Ref. |1|. 

Here we deal with the remaining class of popular stochastic models, in which the 
diffusion anomaly stems from the explicit position or time dependence of the diffusion 
coefficient. Modelling anomalous diffusion via a coordinate dependent diffusivity goes 
back to Richardson’s analysis of relative diffusion in turbulent flows with k = 3 [5]. 
Systematic variations of the local diffusion coefficient in space were, for instance, 
demonstrated in the cytoplasm of living biological cells [22] • A series of publications 
recently explored the stochastic properties of such heterogeneous diffusion processes 
(HDPs) [3ni[3Tl[321 |33l|3ll[35| . Instead of the position dependence, Batchelor introduced 
a time dependent diffusion coefficient for the description of Richardson diffusion [35] . 
Scaled Brownian motion (SBM) with a power law form for the diffusivity is very popular 
in the phenomenological description of anomalous diffusion [37]. Its stochastic properties 
were analysed in detail in Refs. [301I3E1E91I1Q111IIS2]. SBM was shown to grasp essential 
features of granular gases in the homogeneous cooling phase [33] . 

Physically, HDPs and SBM are quite different; HDPs are truly Markovian 
processes (with multiplicative noise [301 [3T] 1 while the time dependence of the diffusion 
coefficient in SBM already heralds its fundamentally non-stationary character. Indeed 
experimentally it is often not possible to separate the effects of time and position 
variation of the diffusivity—compare, for instance, the data in Ref. [H]. This poses 
the question how exactly the position and time dependence of the diffusion coefficient 
built into HDPs and SBM conspire. How do they compete with each other? On the 
level of the stochastic Langevin we here analyse in detail generalised diffusion processes 
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D{x,t) ~ (l-f 

where the physical dimension of the prefactor Dq is [Dq] = cm^“"sec 


^ The choice (|^ 


(2) 


thus combines the two parent processes (HDP and SBM) multiplicatively. In particular, 
we derive the ensemble and time averaged MSDs as well as their statistical properties. 
We also analyse the ageing behaviour and the effects of external conhnement. We will 
discuss similarities and disparities with other anomalous stochastic processes. 

2. Observables 

We start by introducing the physical observables, that we are going to analyse in the 
remainder of this work. The most standard quantity to classify a stochastic process is 
the MSD ([^, which is defined in terms of the spatial average 



(3) 


over the available space region, where P{x,t) is the probability density function to 
find the test particle at position x at time t. This is typically a good measure when 
many short trajectories of tracer particles are available. In many single particle tracking 
experiments, however, few long time traces x{t) of length T (the measurement time) are 
available. The time series x{t) is usually evaluated in terms of the time averaged MSD 



(4) 


where A is the lag time. Often, the additional average 



(5) 


i=l 


over N individual trajectories is taken to produce smooth results. 

For Brownian motion, the increments are stationary, and the time averaged MSD 
in the limit of long times becomes identical to the ensemble MSD; limT/A-^^oo (5^(A) = 


(a;^(A)) in im SHI 07]. Due to this identity we call Brownian motion an ergodic process 


08] . The same is true for anomalous diffusion processes driven by fractional Gaussian 
noise [261119],[50] (although pronounced transients may become relevant [Ill El]), as well 
as for diffusion on fractals [S^. In contrast, non-stationary processes are subject to the 
disparity 


lim (52(A) 7^ (a;^(A)). 

T/A^oo 


( 6 ) 


This behaviour is referred to as weak ergodicity breaking [11051061E71E31ESE51E6]. 
For a range of anomalous diffusion processes, the time averaged MSD in the limit A <^T 
was found to scale linearly in the lag time. 



J'l-K ’ 


(7) 
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despite the anomalous scaling of the ensemble MSD Q. Simultaneously, the explicit 
dependence on the measurement time shows that the process is progressively slowing 
down (0 < K < 1) or picking up speed {k > 1). The specific form ([^ was derived, 
inter aha, for subdiffusive continuous time random walks [5ll ESI [56], HDPs [301 El]) 
and SBM [301130] 130] . We note that a similar linear lag time dependence of the time 
averaged MSD was also observed in ultraslow processes with a logarithmic form of the 
ensemble averaged MSD [121 EH ESI ESI 100] . 

Ergodic processes are reproducible in the sense that each time we evaluate some 
observables for a—sufficiently long—measured time series we obtain the same result 
with only minor deviations. Quantified in terms of the dimensionless variable [30] 





( 8 ) 


this means that the associated distribution 0(^) converges to the delta function, 
lim'T/A^.oo 0(0 = — 1) [H SSI SSI ES]. Weakly non-ergodic processes have different 

limit distributions for 0(0- For instance, subdiffusive continuous time random walks 
have a finite value 0(.^ = 0) and a distinct contribution away from the ergodic value 
^ = 1 [H SSI ESI ES] . HDPs have shapes of 0(0 that are similar to Gamma distributions 
with 0(^ = 0) = 0 [31]. SBMs are weakly non-ergodic but asymptotically reproducible 


The spread of the time averaged MSD at given lag time A ^ T is quantified by 
the ergodicity breaking parameter [H SSI SSI ES] 

EB(A) = ({^(A)) - «(A)>^ = ({=(A)) - 1. (9) 

For Brownian motion EB vanishes linearly with A/T [1], 


4A 

lim EBbm(A) = —. 

A/T-s-O Ol 

Sometimes, the alternative ergodic parameter 


TH(A) = 


52(A) 


( 10 ) 


( 11 ) 


(x2(A)) 

is invoked to provide some additional information about the ergodic properties of the 
process [aisi]. 

The ergodicity breaking parameter ([^ for any finite ratio A/T for non-ergodic 


processes is larger than the corresponding Brownian value (10) and may not vanish in the 


limit of long measurement times. As we show below the ergodicity breaking parameter 
0 is a sensitive measure for the physical processes generating a given anomalous 
diffusion process. Other observables such as the ensemble and time averaged MSDs 
often feature similar or even identical functional forms for different processes and are 
thus not suitable to discern a specific process, compare the discussion in Refs. [31 162] l63]. 

Finally, we address the ageing behaviour of the anomalous diffusion processes. 
This is the explicit dependence of physical observables on the ageing time ta elapsing 
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between the initiation of the process at t = 0 and the start of the measnrement at 
ta- In particnlar, for the time averaged MSD of the aged system the evalnation of the 
associated time series is then performed in terms of 




x(t + A) 


x{t) 


dt. 


( 12 ) 


ta 


For snbdiffnsive continnous time random walk processes [6l], HDPs [3l], and SBM [H], 
it was fonnd that the aged and non-aged time averaged MSDs for A -C T fnlhl the 
relation 

(^2(A;«J)=aAA (13) 

and thus differ only by the ageing factor 



containing, multiplicatively, all information on the ageing and measurement times. 

In what follows, we analyse these quantities in detail for the GDP composed of the 
HDP and SBM parent processes. 


3. Generalised diffusion processes 


Based on the spatio-temporal dependence (|^ of the diffusion coefficient we define the 
generalised diffusion process in terms of the overdamped Langevin equation 

^ = ycpyj X e((). ds) 

The noise is Gaussian with unit variance and zero mean [311 HO] . For the position 
coordinate x{f), this is a multiplicative equation, which we extensively tested for HDPs 
[311 ESI |33l [3l]. With respect to the time dependence of D{x,t), we analysed this 
Langevin equation in the context of SBM processes [ID]. 

We interpret the Langevin equation in the Stratonovich sense. In the mid-point 
discrete version, the particle displacement in the {i + l)st step thus becomes [DT] 


Xi+i -Xi = A 2 




X ( 2 /i+l 2/i), 


(16) 


where the increments {yi+i — yi) of the Wiener process represent a (5-correlated Gaussian 
noise with unit variance. Unit time intervals separate the consecutive iteration steps. 
To avoid a trapping of particles in the regions of zero diffusivity for certain exponents 
for a and /3, the form (16) chosen for the simulations was regularised by addition of the 
small constant Dos = 10“^, compare Refs. 


for more details on the simulation 
procedure. The particle’s initial starting position in all the results presented below is 
Xq = 0.1 (see Ref. [M] for details of the effects of Xq). The discrete scheme (16) is 
supplied with either natural or reflective boundary conditions in the following. 
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From the Langevin equation (15) it is straightforward to derive the diffusion 
equation (compare Ref. [3T] ) 


d 


d 




D(x,t)P{x,t) 


(17) 


Following the Stratonovich interpretation we see that the diffusion coefficient appears 
symmetric with respect to the Laplacian operator. In the limits /3 —)■ 0 and a —)■ 0, 
respectively, we recover the diffusion equations with purely position dependent [31] and 
time dependent HD] diffusivity. 

3.1. Unconfined motion 

For unconfined motion we apply the natural boundary conditions hm| 2 ;|^oo P{x, t) = 0 

ohc 

2; 12/p 


to the GDP diffusion equation (17) to find the probability density function 

|^|i/p-i 


P{xfi) = 


^/A^TDfi^+d 

where we used the abbreviation 
2 

P = 


exp 




2 — a 


(18) 


(19) 


The functional form of the probability density function (18) thus combines the spatial 
features of the HDP process with the space dependence ~ |x|" of the diffusivity of SBM 


and the modified time dependence due to the contribution ~ of the diffusivity 
HD]. In particular, for p > 1 (p < 1) we observe a stretched (compressed) Gaussian 
shape for larger x as well as a cusp (dip to zero) close to the origin. This behaviour is 
supported from simulations in Fig. 


Integration based on Eq. (18) produces the ensemble averaged MSD 
{xfit)) = 

The scaling exponent 


F(l/2 + p) 


TT 


1/2 


2p 

1) 


( 20 ) 


k = (1 + /3)p (21) 

is thus given by the product of the HDP and SBM exponents. In the limits (3 = 0 and 
a = 0 we recover the scaling exponents n = p and k = 1 + (3 oi the pure parent processes 
for HDP and SBM, respectively [ani EH HD]. For the parameter range —2 < a <2 and 
—1 < /3 < 1 of the HDP and SBM processes, the phase space of the GDP is depicted 
in Fig. 1^ for all parameter values in the red (blue) areas, the GDP is super diffusive 
(subdiffusive). Note that we also simulated the GDP for the temporal exponent (3 = ?>/2 
outside of the typical SBM range — 1 < /3 < 1 in Fig. 1^, observing good agreement 


with the prediction for the MSDs (20) and (22). 


In analogy to the derivations in Refs. [miiQ], we obtain the time averaged MSD 
in the limit A <tlT, 

F(l/2+p) /2\^^ „„ A 


52(A)\ 


TT 


1/2 


2\ 2P 

pj 
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Figure 1. Probability density function P{x,t) for free and confined (see below) GDPs 
for (a) subdiffusive and (b) superdiffusive choices of the critical exponents a and /3, as 
indicated in the plots. For subdiffusive GDPs the small spike at x = 0.1 is a remainder 
of the initial starting position of the particles. For superdiffusive GDPs this is not 
visible due to the fast formation of the peak of P{x,t) at x = 0 in the region of the 
slowest diffusivity. Parameters: T = 10^, N = 10^. The L values are as indicated. 



Figure 2. Phase space of the generalised diffusion process. In the red (blue) areas 
the GDP is superdiffusive (subdiffusive). The dashed lines at a = 0 (/3 = 0) represent 
pure SBM (HDP), respectively. 
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Figure 3. a) Ensemble and time averaged MSDs for unconfined GDPs for a = —1/2, 
P = 3/2, T = 10®, and N = 100. b) Ensemble and time averaged MSDs of unconfined 
GDPs with a = 1 and /3 = —1/2, that is, for a linear time dependence of the ensemble 
MSD (k = 1). Note the initial relaxation of the ensemble MSD to the asymptotic 
behaviour {x^{t)) ^ E. The inset shows the behaviour of the ergodicity breaking 
parameter with the Brownian asymptote (10). c) Ensemble and time averaged MSDs 
for unconfined SBM process for /3 = 3/2. The limiting laws (20) and (22) with p = 1 
are represented by the dashed curves. Parameters: T — 10®, N = 50. 


hhA)) (^ 


l-(l+^)p 


( 22 ) 


Time and ensemble averages are thns disparate, and in analogy to the parent processes 
HDP and SBM we obtain weak ergodicity breaking. The alternative ergodicity breaking 
parameter for the GDP thns becomes 




(23) 


The explicit dependence on the measnrement time T in relations (22) and (23) 
is a signatnre of the non-stationary character of the GDP. As fnnction of T the 
amplitnde of the time averaged MSD continnonsly decreases (increases) for snbdiffnsive 
(snperdiffnsive) parameters. 

Fig-i shows the time averaged MSD from simnlations for a range of valnes of the 
scaling exponents a and fd. The predicted scaling (S^) ~ according to Eq. (22) 
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Figure 4. Time averaged MSD for GDPs in the limit A/T ^ 1 (A = 1). The 
asymptotes are given by Eq. (22) and the values of the exponents a and /3 are as 
indicated. Good agreement with the predicted scaling behaviour is observed. 


agrees nicely with the simulated data. The linear lag time dependence (5^) ~ A is 
supported by the simulations shown in Fig. and below in Fig. 

Note that for those combinations of the scaling exponents a and (3 for which the 
ensemble MSD grows linearly [k = 1, see the diagonal separating subdiffusion and 
superdiffusion in Fig. the GDP is still weakly non-ergodic, as illustrated by the 
irreproducibility of the individual time averaged MSDs and their inequivalence with the 
ensemble MSD shown in Fig. for a = 1 and f = —1/2. Despite the fact that the 
trends of the spatial and temporal variation of the GDP diffusivity compensate each 
other in terms of the MSD the overall process is still non-stationary. 

The scatter of the time averaged MSD from individual simulated trajectories is a 
sensitive function of the diffusivity D[x, t), as shown in Fig. The spread of S'^ becomes 
broader as the value of the spatial exponent a approaches its critical value a —?• 2. This 
is a characteristic property of HDPs |31]- With increasing A the spread of traces 
increases only slightly. For subdiffusive GDPs the distribution 0(,^) is similar to the 
asymmetric Rayleigh form discussed in Ref. [HT] , 

The ergodicity breaking parameter EB(A) of GDPs is illustrated in Fig. [^. The 
deviations from the Brownian result ( |Io| ) in the limit A/T 1 quantihes the departure 
of the diffusing particles from the ergodic behaviour. We hnd that for GDPs the value 
of the ergodicity breaking parameter does not vanish in the limit A/T -C 1, see Fig. &■ 
For GDPs the value of EB in the limit A/T —0 varies signihcantly with /3, as shown 
in Fig. [^. We note that the diffusing particles under lesser conhnement produce larger 
values of EB, compare the curves in Fig. |^, see below. We also note that varying initial 
conditions xo alter the values of EB (not shown). The auxiliary ergodicity breaking 
parameter for free and conhned GDP motion is shown in Fig. |^. 
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Figure 5. Amplitude scatter distribution (j>{E} of individual time averaged MSD traces 
of the GDP for varying lag times A and the extent of confinement. As parameters we 
chose T= 10^, N = 10^, as well as for the top row a = —2 and /3 = 1/2 (subdiffusion), 
for the bottom row a = l,/? = l/2 (superdiffusion). Different degrees of confinement 
were imposed, the left panel in each row corresponds to effectively free motion. 


3.2. Ageing effects 


The effect of ageing, that is, the dependence of physical observables on the time 
difference ta between system initiation at t = 0 and start of the measnrement of the 
process in analogy to the results for HDPs [M] and SBM SB factorises according to 


Eq. (13). The ageing factor (14) 


I ^ 


, N p(l+/3) 




(24) 


now has a parametric dependence on both p = 2/(2 — a) and fd. This scaling with 
the ageing time is indeed supported by our computer simulations for a number of 
combinations of the scaling exponents a and //), as shown in Fig. 0 Fig.0, shows 
that for aged GDPs the scaling function (24) remains valid also for // > 1, see the 
discussion for ageing SBM below. Note that for subdiffusive GDPs in general longer 


traces are needed to obtain the same degree of convergence to Eq. (24), as detailed in 

Fig. [ 5 . 
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Figure 6. Ergodicity breaking parameters EB (panel a) and 8B (panel b) for free 
(L = 10®) and confined superdiffusive GDPs (/3 = 1, a = 1/2). The asymptotes (10) 
for Brownian motion in (a) and (23) in (b) are shown as the dashed lines. The values 
of the confinement L are as indicated. Parameters: T = 10®, N = 150. Panel c): log- 
linear plot of EB at A = 1 for different /3 values for free GDPs. Parameters: T = 10"^, 
N = 500. 


Ageing SBM. Let us briefly digress to study the ageing behaviour of pure SBM. From 
the results of Ref. BDl together with the definition of the ageing time averaged MSD 
(12) we find the full form of the ageing factor 


mJT) = 


(1 + + {tJT)f‘+'‘ - (tJT + A/r)f+2 - (1 + tJT - 

i-(A/r)'>+2-(i-A/ry+2 


(25) 


that for short ageing ta/T <C 1 and lag times A/T <C 1 reduces to relation (24) with 
p = 1. This implies that for superdiffusive SBM with /I > 0 the magnitude of time 
averaged MSD increases with the ageing time ta/T and decreases for subdiffusive SBM 
with /3 < 0, as shown in Fig. Physically this is a consequence of the continued 
acceleration of particles in the superdiffusive case, contrasting the progressively localised 
particles for subdiffusion. We refer to Ref. HU for more details on the properties of 
ageing SBM including the mean first passage time density. 

We support the limiting form (24) with p = 1 of the ageing factor for ageing 
SBM from computer simulations in Fig. for a range of scaling exponents The 
agreement of the simulations with the asymptote (24) is particularly good for the case 
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Figure 7. Ageing factor Kp^pita/T) for aged GDPs in the limit A/T <C 1 (A = 1). 
The values of a are given in the plots. We show the values /? = 5/6 (panel a), /? = 3/2 
(panel b), and /3 = —1/2 (panel c). The asymptotes (24) are represented by the dashed 
lines. Parameters: the overall length of the traces is 10^ steps, of which the evaluated 
trace lengths T are indicated in the plots. We averaged over N = 200 traces for each 
data point. 


of superdiffusive SBM with T = 10^ steps. Strongly subdiffusive SBMs likely require 
longer traces to converge better to (24). Performing simulations for the range /? > 1 we 
also reveal excellent agreement with the ageing factor (24), see the top curves in Fig. 


3.3. Confined motion 

Let us now address GDPs in the conhnement of hard walls. To this end we establish 
reflecting boundary conditions at a; = ±L. Fig. demonstrates that for weak 
conhnement (large L values) and at shorter (lag) times the ensemble and time averaged 
MSDs behave similarly to those of the unconhned process, as expected. At stronger 
conhnement and at longer (lag) times the MSD saturates at the plateau value 


( 26 ) 
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Figure 8. Ageing factor A/j for ageing SBM, the dashed lines represent the asymptote 
(24) for p = 1. Parameters: the total trajectory length is 10® steps and the length T 


evaluated in the graphs is stated in the plot panels. The averaging was over N « 100 
trajectories for each /3 value, and we chose A = 1. 


The ct-dependence of the value (26) in the prefactor p~^/^ reflects the fact that the 
associated stationary probability density function is still position dependent due to 
the variation of the diffusivity, see also Fig. For a = 0, the process assumes 
an equidistribution with the value 1/(2L). We stress that the plateau value (26) is 
independent of the temporal scaling exponent f. In the hard conhnement by inhnitely 
steep walls the stationary state is generally independent of any ^-independent noise 
strength. The variance (26) is thus geometry induced and not due to the competition 
of thermal noise and confinement as, for instance, in an harmonic confinement. In the 
latter case, no stationary solution is reached for SBM m. 

Fig. also shows that for the same interval width L subdiffusive GDPs require 
longer times to reach the stationary values. After reaching the stationary plateau, 
the value of the time averaged MSD is twice that of the ensemble MSD, an effect 
of the very definition of the time averaged MSD |Tj. Note that due to the pole at 
A —)■ T in the definition (|^ of the time averaged MSD the asymptotic equivalence 
limA-s.T(<5^(2^)) —t {x^{T)) holds |1]. 

The probability density function P{x, t) of confined GDPs for different values of the 
scaling exponents a and f3 as well as of varying degrees of confinement is shown in Fig. [Tj 
It is seen that when the particle is reflected repeatedly in stronger confinement by the 
walls, the tails of P{x, t) become perceivably raised, to fulhl the boundary condition 
dP{xp)/dx\^^f^ = 0 . 

The distribution of amplitudes of individual time averaged MSDs described by 
0(^) for confined GDPs is shown in Fig. The width of 0(^) decreases for more severe 
conhnement, as expected. The width of 0(^) for shorter measurement times T is larger 
(not shown), due the smaller number of rehections from the walls. For conhned GDPs 
becomes more localised and symmetric both for subdiffusive and super diffusive 
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Figure 9. Ensemble and time averaged MSDs for confined subdiffusive (a = —2 and 
P = 1/2, top row) and super diffusive (^0 = 1 and a = 1/2, bottom row) GDPs. The 
limiting behaviours (201 and (22) for unconfined motion as well as the plateau value 
(26) under confinement are indicated by the dashed curves. Parameters: T = 10® and 
N = 150. The width of the confining interval L is indicated in the panels, the value 
L = 10® corresponds to an effectively unconfined situation. 


combinations of the scaling exponents a and /3, as demonstrated in Fig. |^: as one may 
already venture from Fig. the spread of individual 6^ decreases severely when extreme 
excursions of particles causing large trajectory-to-trajectory variations are progressively 
impeded by the conhnement. 

The dependence of the ergodicity breaking parameter on the lag time A for different 
interval lengths L is shown in Fig. [^, and the variation of EB with the measurement 
time T is illustrated in Fig. We hnd for GDPs the reciprocal dependence 

EB(T) ~ t. (27) 


with the trace length T, similar to the behaviour of conhned HDPs [Mj . Similar to HDPs 
and SBM, for free and conhned GDPs our simulations conhrm Eq. (23) in Fig. [^. In 
the long time limit the saturation plateau is reached and £B{A) —2. We point out 
that this behaviour is signihcantly different from that of subdiffusive continuous time 
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Figure 10. Ergodicity breaking parameter EB for confined GDPs from simulations. 
The asymptote (271 is indicated by the black dashed curve. Notations for the symbols 
are the same as in Fig. 


random walks, for which the power-law ((5^(A)) ~ [A/Tf " is obtained |53j . 

4. Discussion and Outlook 

We studied the non-homogeneous and non-st at ionary generalised diffusion process with 
the power-law diffusivity D{x,t) ~ We explored analytically and by extensive 

computer simulations the competition between the two parent processes revealing a 
number of universal characteristics. While the general scaling forms for the ensemble 
and time averaged MSDs are similar to those of the parent processes HDP and SBM— 
and in fact to those of subdiffusive continuous time random walks in im EH]— we 
showed that for those combinations of the scaling exponents a and fd which produce 
a unity anomalous diffusion exponent, k = 1, the non-st at ionary character of the 
GDP is still visible in the irreproducibility of the time averaged MSD. For the ageing 
GDP we obtained that the ageing factor has the same functional form as for HDP, 
SBM, and subdiffusive continuous time random walks miMisi]. We also explored the 
properties of GDPs conhned in an interval including the saturation of the MSD and the 
variation of the ergodicity breaking parameter with the trace length and the degree of 
conhnement. Physically, the combination of temporal and spatial variations of the local 
diffusivity appears both a natural and attractive concept to model anomalous diffusion, 
for instance, in microscopic biological systems. 

Restricted diffusion in various porous media was considered in terms of effective 
time dependent diffusivities. Such approaches may provide some information about 
the typical size of restricted regions or cavities experienced by a tracer. The concept 
of time dependent diffusivities is, for instance, applied to interpret signals in NMR 
experiments [651166| l67j aimed at determining the degree of tissue interconnectivity and 
permeability. Water and ion diffusion in various human tissues is often anisotropic and 
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sometimes anomalous [HE]- Anisotropies in water diffusion were observed in biological 
tissues with a directional/fibrous structure, for instance, in the white matter of the 
human brain, nerve fibres, and muscle fibres [66]. 

Specifically, the instantaneous diffusivity D{t) = d{r^(t))/{6dt) in muscle fibres was 
shown to assume the long-time limit to decay D(t) = D{oo) + const/t^A_ Alternative 
power-law scaling of the form D{t) ~ corresponding to different structures of the 
disorder in the system were also examined, p often ranging in the interval 1/2 < /r < 3/2 
[66j. The value of p = 3/2 is known for a fully uncorrelated medium in three dimensions. 
In a biological context, the power /i of the tail of the diffusivity depends on the 
positioning of boundaries in the system. The latter are typically cell membranes that are 
only weakly permeable to water or ions. It was argued that some deceases may caused 
changes in the brain or muscle tissues that can be detected by measuring the changes 
in the time dependent diffusivity and the exponent p as compared to healthy tissues. 
Note that in these biological examples with time-dependent diffusivity no complete 
localisation of particles in the long-time limit exists and the basal diffusivity D{oo) is 
hnite. 

In porous interconnected media, for instance, in a suspension of reflecting spheres 
[69j, the long-time diffusivity of a tracer can be defined as D(t) = (r^(t))/(6t). It 
can again be represented in the time dependent form D{f)/DQ Ai + Aj/t + Ag/e/'^ 
1701 . The leading time dependence for the instantaneous diffusion coefficient is then 
D{t)/DQ ~ Ai —O.SAst Here 1/Ai is the tortuosity of the medium that renormalises 
the long-time tracer diffusivity due to the presence of inter-connected cavities |69j . 

Are such assumptions of an effective, time dependent diffusivity always justified? 
In fact, any anomalous diffusion law ([^ can be interpreted in terms of the effective 
diffusivity A'eff(t) — such that we can rewrite Eq. Q in the form {x^{t)) ~ Kes{t)t. 
However, this is not sufficient to identify the observation of anomalous diffusion with 
the GDP or SBM models studied here. It is important to either have more detailed 
knowledge of the underlying physical process, which may support the GDP approach 
or rather other physical models such as continuous time random walks [U 06]. In 
particular, we note that despite many close similarities, SBM and fractional Brownian 
motion (FBM) capture physically strongly different situations. If direct physical insight 
is not possible, the stochastic properties of the process need to be scrutinised in more 
detail by using complementary diagnostic tools 1111621 in]. 

Similar to these examples, position dependent diffusivities were explored in a 
range of systems including tracer diffusion in groundwater aquifers [7|. The position 
dependence is directly accessible experimentally and was systematically sampled for the 
motion of smaller protein tracers in living biological cells [29] . 

We note that HDP and SBM processes may also give rise to ultraslow diffusion with 
a logarithmic dependence of the MSD. For the HDP this corresponds to an exponential 
variation of the diffusivity in space [32] while for SBM ultraslow motion emerges in the 
limit /? —)■ — 1 [12]. Such ultraslow diffusion is equally weakly non-ergodic and ageing, 
in analogy to diffusion in ageing environments [58], diffusion in quenched Sinai disorder 
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